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A NEUTRAL PARTICLE WAKE EXPERIMENT FOR MEASURING UPPER 
ATMOSPHERE TEMPERATURES FROM A MOVING VEHICLE 

L. H. Brace and W. C. Wu 
ABSTRACT 

A method is discussed for measuring the kinetic temperature of the neutral 
particles in the upper atmosphere. A suitable baffle is mounted an appropriate 
distance in front of the orifice of a satellite-borne neutral particle mass spec- 
trometer. As the satellite spins, the neutral particle wake formed by the baffle 
scans across the orifice and causes a bite-out of the pressure maximum which 
would otherwise be observed by the spectrometer as the orifice rotates through 
the velocity vector. At very low temperatures the baffle would completely block 
the orifice from the stream of thermal particles through which the satellite is 
moving, and the pressure in the spectrometer would be reduced to zero as the 
orifice passed through zero angle of attack. At atmospheric temperatures (500- 
2000°K), a suitable baffle produces a pressure reduction which varies approxi- 
mately inversely with the temperature. 

Since the various neutral particle species have different thermal velocities, 
each will independently produce its own wake. The baffle performs a kind of 
mass discrimination in which the lateral components of thermal motion cause the 
lighter species to form the near-wake and the heavier species to form a more 
extended wake. Thus, at any given distance behind the baffle, a greater pressure 
reduction is observed for the heavier species. 

When chemically active species such as atomic oxygen are sampled the depth 
of modulation will still reflect the temperature, but the absolute amplitudes of 
the pressure inside the chamber will depend upon the amount of recombination. 
When mixtures of atomic and molecular species are present in the atmosphere, 
this technique provides an additional source of information regarding the relative 
concentrations of the atomic and molecular species. 
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A NEUTRAL PARTICLE WAKE EXPERIMENT FOR MEASURING UPPER 
ATMOSPHERE TEMPERATURES FROM A MOVING VEHICLE 


INTRODUCTION 

The gas kinetic temperature (T Q ) is a parameter of considerable interest in 
the physics of the upper atmosphere. Nearly all of our knowledge about its be- 
havior has been inferred from neutral particle scale heights based on satellite 
drag studies (Jacchia, 1964). Direct measurements from satellites have not 
generally lead to entirely reliable neutral particle scale heights because of 
latitudinal and local time structure which is superposed upon the changes of 
neutral density with altitude. Therefore, reliable measurements of gas temper- 
ature have been limited to a few rocket flights in which true vertical profiles of 
individual species have been measured (Spencer et al, 1965) (Nier, et ai, 1964). 

The major difficulty encountered in measuring T directly from i. satellite, 
or a rocket, is that the translational velocity of the measuring platform greatly 
exceeds the thermal velocity of the particles. Therefore the flux of particles 
flowing to a sensor experiences extremely wide variation with angle of attack. 

At satellite velocities, these variations are relatively insensitive to T Q . The 
method discussed here overcomes this difficulty by interrupting the directed 
stream by means of a rectangular baffl j. The method is very sensitive to T 0 
because the particles which populate the wake of the baffle do so only because 
of their thermal motion. It will be shown that the amplitude of the pressure in 
the wake, or conversely the depth of pressure reduction in the wake, is a strong 
function of the kinetic temperature. 


DESCRIPTION OF THE EXPERIMENT 

If the measurement is to be performed on a spinning satellite, as we have 
thus far assumed, a long rectangular baffle is mounted at a distance of several 
inches in front of a small orifice leading to the ionization chamber of a neutral 
particle mass spectrometer. This arrangement is mounted perpendicular to the 
spin axis of the spacecraft as shown in figure 1. As the satellite spins, the 
neutral particle wake which forms behind the baffle scans across the orifice 
and reduces the flux of neutrals into the chamber. This causes a reduction of 
the partial pressure of each species as the angle of attack passes through zero. 
At the limit of zero temperature, the baffle would completely block the orifice 
from the stream of particles through which the satellite was moving, and no 
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BAFFLE 



Figure 1-Experimental arrangement. As the satellite spins, 
the long rectangular baffle (seen from the narrow dimension) 
passes between the velocity vector, w, and the orifice. The 
wake which forms behind the baffle scans across the orifice 
and causes a pressure reduction which depends upon the ther- 
mal velocity of the particles sampled. 

thermal particles could enter the orifice. However, at the high temperatures 
found in the upper atmosphere, the reduction in pressure behind the baffle is 
much less because the thermal motion of the particles causes them to fill-in be- 
hind the baffle more rapidly. This effect is illustrated in figure 2 which shows 
the relative sizes of the wake for H , H e , 0 and N 2 behind a very large baffle. In 
constructing this drawing we have assumed that the thermal velocity of H is 
comparable to the satellite velocity. Clearly the depth of the pressure well in 
the wake increases with the mass of the particle. The depths of the "pressure 
wells "for a fixed translational velocity depend upon the choice of R/ 1 , the prime 
variable in designing the experiment. 
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SPACECRAFT SURFACE 



Figure 2-Relative sizes of the wake for various neutral 
species, assuming the velocity of the satellite is com- 
parable to the thermal velocity of atomic hydrogen. The 
wake boundaries are arbitrarily defined. 
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ANGLE OF ATTACK ( 8 0 ) 

Figure 3-Chamber Pressure of N 2 as a function of angle of attack, for R/'i = 0, 0.05 and 0.1. 



CHAMBER PRESSURE (ARBITRARY UNITS) 



Figure 4-The effect of atmospheric temperature upon the pressure well 
for a baffle R/^ = 0.U5. The temperature affects both the depth and 
width of the pressure trough. 
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Figure 5-The effect of atmospheric temperature upon the pressure well for a baffle R/'i = 0.1, 
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scale is expanded about zero to show only the portion of the pressure character- 
istic which is affected by the baffle. Nearly the entire range of exospheric 
temperature is encompassed by the three curves shown (500°, 1000°, 2000°). 

The atmospheric temperature modifies both the depth and the width of the 
pressure well. The pressure reduction ’’atio P . /P , is 0.60 for 2000°K and 

m in max 

0.265 for 500 K when the 0.05 baffle is used. When the 0.1 baffle is used, the 
pressure reduction ratio is 0.25 for 2000°K and 0.022 for 500°K. From this it 
can be seen that, in theory, the larger baffles provide greater resolution in the 
temperature. In practice, however, the minimum pressure which can be resolved 
depends upon experimental factors such as the background pressure due to con- 
taminants, the gas time constant of the orifice-chamber system, the sensitivity 
of the spectrometer and the telemetry resolution. 

Although the background pressure can be measured during each spin cycle 
as the orifice looks back into the spacecraft wake, the background itself may be 
orientation dependent, or time dependent. 

The gas time constant of the system can be made adequate by maintaining a 
large ratio of orifice area to chamjer volume as well as by selecting a suitably 
low spin rate. Spin rates of as great as 20 rpm should be possible. The effects 
of the orifice radius upon the pressure characteristics are second order, as 
discussed in Appendix B. Therefore the gas time constant can be accommodated. 

The sensitivity of the spectrometer ultimately limits the altitude to which a 
given species can be employed for the neutral wake measurements. The most 
abundant species are the most desirable in this regard, therefore N 2 and 0 are 
the prime constituents to be employed in the Earth's atmosphere. 

The telemetry systems currently in use provide telemetry accuracy of about 
1/2% of full scale. Although current detection schemes could be devised for in- 
creased electronic resolution of the pressure well (P min ) * the telemetry accuracy 
will ultimately limit the depth of the pressure well which could be resolved. 

For the reasons given above the R/'t should be selected to produce perhaps 
less than a factor of ten reduction in the partial pressure of the species of 
prime interest. 


The Pressure Curves for Other Particle Masses 


Although all of the calculations represented in figures 3, 4 and 5 were 
carried out for N 2 , some feeling for the effect of employing other species can 
be obtained from these figures because the thermal velocity ( c m ) depends upon 
both the temperature and mass of the particle. 
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c 

m 


( 2kT 0 /m) 1 2 


Thus a factor Oi our increase in T 0 is equivalent to a factor of four decrease 
in m . For example, in a theoretical sense, the curves in figures 4 and 5 labelled 
500°, 1000° and 2000°K could equally well have been labelled masses 28, 14 and 
7, at a fixed temperature of 500°K. 


The Measurement of Atomic Species 


Atomic constituents such as 0 can also be employed for determining the 
temperature. At most satellite altitudes ( >250 km), nearly all of the oxygen 
exists in atomic form. Although the influx of 0 into the orifice will be described 
properly by the theory, the flow of 0 out of the orifice (equation 8) will be re- 
duced by the amount of recombination occurring on the walls of the chamber. As 
a result, the amplitude of the partial pressure curves for 0 will be less than 
would otherwise be expected. However, the relative shape of the O curves is not 
necessarily changed by surface recombination unless the percentage recombina- 
tion is pres sure -dependent. It seems reasonable to assume that the recombina- 
tion factor would be constant as long as the mean free path is much greater than 
the chamber dimensions. In this case, the flux of 0 into the chamber (G. ) is 
balanced by the out flux (G 0 ) and the recombination flux (oXT). 


or 


G i = G 0 + aG i 

G 0 = G.(1 - a) 


where a is the fraction of incoming particles which recombine before leaving. 

If a is constant as we have assumed, then the recombination flux will simply 
change the ratio of the incoming and outgoing 0 flux. This will reduce the ab- 
solute amplitude of 0 pressure but will not change its relative pressure modu- 
lation. Since the temperature is derived from the latter, it should be feasible 
to employ the atomic oxygen data for the measurement of temperature. This is 
important because atomic oxygen is the major constituent at the altitudes where 
this technique would be applied. 

From the above discussion it is apparent that any attempt to determine the 
absolute concentration of atomic constituents requires that the molecular form 
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also be measured so that the recombination losses can be evaluated. If the only 
source of 0 2 is the recombined 0 , then both the 0 and 0 2 pressure character- 
istics carry the same temperature information. Thus the temperature can be 
obtained from either without prior evaluation of the recombination fraction, if 
this fraction is independent of the chamber pressure. 


Mixtures of Atomic and Molecular Species 

If the neutral wake method were employed at lower altitudes where both 0 
and 0 2 are abundant, (as would be encountered below 200 kilometers on rockets 
or low perigee satellites) the arguments of the previous paragraphs could be 
inverted to determine the relative concentrations of atomic and molecular oxygen 
in the atmospheric sample. In this case, the temperature would be derived from 
the N 2 (or 0) pressure characteristics. The pressure modulation in 0 would 
still reflect only the gas temperature, but the 0 2 pressure characteristics would 
now have two components (1) the internally recombined oxygen having the same 
shape as the 0 characteristics, and (2) the atmospheric 0 2 having a much greater 
depth of modulation corresponding to its greater mass. Thus the relative concen- 
trations of 0 and 0 2 in the atmosphere control the relative depth of the pressure 
well observed in 0 2 . In practice the relative concentrations of 0 and 0 2 would 
be derived by adjusting their values to reproduce the observed 0 2 pressure well 
theoretically. 


Tne Effects of Particle Baekscattering 


A potential source of error in the neutral wake method arises from the 
double scattering of particles from the satellite surface to the back of the baffle 
and then into the orifice. This effect must be rather small at small angles of 
attack for the following reasons. The same flux of particles which is occulted 
by the front of the baffle also impacts the entire forward surface of the space- 
craft. An equivalent flux is therefore flowing away from the satellite, and some 
of this impacts the inward side of the baffle. Fortunately little of the flux which 
is re-emitted by the baffle actually reaches the orifice because of the small 
solid angle it subtends when viewed from the baffle. In any event, the random 
flux scattered into the orifice will be extremely small compared to the directed 
stream which the baffle intercepts. 

If a given application requires further reduction of the backseattered com- 
ponent, one could shape the inside surface of the baffle to direct the scattered 
particles away from the orifice. The entire experiment might also be mounted 
on a boom to reduce the amount of surrounding area for the initial baekscattering. 
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Some effects of backs catte ring should be evident at angles of attack beyond 
90° where the orifice is looking into the spacecraft wake but the inside of the 
baffle is still receiving the full stream, as shown in figure 6. If these fluxes 
produce detectable pressure, they may be valuable for estimating the amount of 
backscattering occurring at zero angle of attack. 


* 

/ 



Figure 6-A situation in which back scattered particles could 
be observed in the absence of the direct stream. 




In the case of atomic species, the back of the baffle and the surrounding 
area might be treated to enhance the surface recombination, thus minimizing the 
amount of the atomic form which would reach the orifice through backscattering. 


SUMMARY 

In summary, it appears feasible to employ the wake behind a moving baffle 
to investigate the thermal velocities of various neutral particle constituents of 
the upper atmosphere. If the partial pressures of the various species are re- 
solved by a neutral particle mass spectrometer, the temperature of the atmos- 
phere can be derived independently from more than one species. When both 
atomic and molecular forms of a gaseous element are present in the region, the 
depth of the wake provides additional information on their relative abundances. 

We have demonstrated the effects of the following factors: (1) Baffle size 
and mounting distances, (2) atmospheric temperature and mass, and (3) orifice 
size (Appendix B). The effect of particle scattering has been estimated to be 
small. 
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APPENDIX A 

Theory of the Neutral Particle Wake Experiment 

List of Symbols 


T Q = gas kinetic temperature of the atmosphere outside 
T = gas kinetic temperature in the chamber 
k = Boltzmann constant 

n = gas concentration in the atmosphere outside 
n . = gas concentration in the chamber 
rr. = mass of the gas particle 
P 0 = gas pressure in the atmosphere outside 
P. = gas pressure in the chamber 
v = thermal velocity of the gas particle 
w = translational velocity of the instrument 
R = half width of the baffle 
i - distance between the orifice and the baffle 
f s = velocity distribution function of particles 

f = velocity distribution function of particles relative to a moving sensor 
r = radius of the orifice 


In this appendix, we outline the procedure for calculating the pressure modu- 
lation caused by a baffle which is mounted in front of an orificed chamber on a 
spinning satellite, as described in the preceding paper. First we develop the 
integral equations describing the flux of particles into and out of the orifice with- 
out a baffle. By equating these flows and employing the gas law, we arrive at an 
integral equation describing the ratio of internal pressure to atmospheric pres- 
sure, P./P 0 . We outline the method of numerical solution of this equation. We 
then introduce the baffle by suitably modifying one of the integration limits on the 
angle of attack integral. 
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The Theoretical Development 


The neutral particles in the atmosphere are assumed to obey the Maxwell- 
Boltzmann velocity distribution, 


f ,( v )= n o 


m 


,3/2 - 


2kT, 


(v ) 2 


2 7rk T„ 


( 1 ) 


The velocity distribution relative to a sensor system moving with velocity w 
becomes 


f(v) =n 0 


m 

277k T„ 


m — . 2 

3/2 


( 2 ) 


We select our coordinate system in such a way that the z- axis is the orifice 
normal and the velocity vector, w, lies in the x-z plane, as shown in figure 7. 
Using vector notation, 


w = (w sin 0 Q , 0 , w cos d Q ) . 


( 3 ) 


The flux density of atmospheric particles reaching the orifice (the origin) 
can be expressed 


-27T fk-n/2 - 

G * i L L n °fe) 


_ m ( V + w ) 

3/2 2kT rt 


v cos^ dqpd(-cos^)v 2 dv • (4) 


The range of thermal velocities which can flow into the orifice can be expressed 


v = (-v sin 0 cos qp, - v sin 6 sin cp, - v cos 6) . (5) 

After the particles enter the chamber, they become thermalized by multiple 
collisions. The flux density of particles which flow out of the orifice is then 
given by 
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Figure 7-The coordinate system employed in the theory. The 
velocity vector, w, always remains in the x*z plane and the 
angle of attack, 6 q , is the angle between the orifice normal 
(z-axis) and w. 



/ 


d 3 v f s (v) v cos 0 


which can be shown to be 


G 


o 



m v 
tk T £ 


vcos<9d<pd(-cos0) v 2 dv 


1/2 


= n. k T. /(27rmkT. ) 



Assuming no source or sink of particles within the chamber, the in flow and 
the outflow will be identical 


G 


o 



(9) 


By invoking the kinetic gas law 


P. 

1 


= n. k T. 

i i 


( 10 ) 


P„ = n n k T„ 


and combining (4), (5) and (7), the pressure ratio between the chamber and the 
atmosphere is found to be 


Pi / T. 

¥ o VT 0 


where 


k l/2 „ /•27T pir / 2 

’LL [ (s ‘ +i)e '" i/2 


q) - 0 8 = 0 


+ 77 1/2 s 1 (l-5 + s 1 )(l+erf Sj) 
sin#cos#d# dcp 


-( w /c m ) 2 s ? 
e e 1 


Sj = — [ sin <9 q sin# cosqp + cos 0 O cos #] 


( 11 ) 


( 12 ) 


c = (2kT n /m) 1/2 , the most probable velocity 


(13) 


Equation (11) is equivalent to a simpler formula 


1 — 


t. y/2 2 

— — 1 f e_s + S77 1/2 (1 + erf s)] 

T o / 


where 


(14) 


w n 

s = — cos 6. 
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developed by Schultz, et al, 1948. However, the integral form (11) permits geo- 
metrical structure which may be in view of the orifice to be considered. The 
difference (11) and (14) arises from the technique of integration. 


The Numerical Integration Method 


The numerical integration is carried out by computer in the following 
manner. Writing (11) in the form 


where 

F(0, 9) 



(15) 


(16) 


6 is divided into 180 intervals 


( r 7 t 2tt rr 

360 360 2 


and the resulting 181 points are named 


(1, 2, 3, 181). 


<p is divided into 360 intervals 



and the resulting points are named 


(1, 2, 3, .... 361) ■ 


18 


The numerical integration is carried out as indicated 


Letting 


P. /T.\ I/2 iso 

"L E lF(e --^ )+F(e ^ v)] ito^- (17) 

^ ~ i = 1 


360 

8(f) = 2^ tF<e,,f) + F(0 1+1 , <p)]JL, 

i = 1 


T y/2 360 


hy ~ H [g( v + 

' ' i = 1 


7 T 


360 


(18) 


= H ( T i’ T 0’ W ’ C m ’ *0> 


(19) 


In test cases, the numerical calculations by computer agrees with the earlier 
form of the pressure equation (14) within 0.1%. 


Introducing the Baffle 

If the baffle takes the form of a long rectangular strip, as shown in figure 8, 
the integration limit on 6 in equation 11 must be modified. The trigonometric 
relationship of the angles shown in the figure are 


Tan 9 

C 



Cos qp = JL 


Q 


(20) 



Tan -1 


Q 

i 


- Tan" 1 


R 

£ 


1 

cos cp 


sin 


R 

[R 2 + cos qp) 2 ] 1/2 
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Figure 8-The baffle coordinate system. 


where 0 is the critical angle below which a particle arriving from azimuth 
angle qp will have no chance of reaching the origin (the orifice). The selection 
of baffle half width, R , and distance, t , thus controls this critical angle which 
becomes the lower integration limit on 6 , as shown in equation (21), 





F ( 0 , <p) dt? dcp • 


( 21 ) 


The numerical integration process described earlier was employed to solve | 

(21) for the several baffle dimensions and atmospheric temperatures represented 
by figures 3, 4 and 5. 

i 

I 
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APPENDIX B 

Consideration of Orifice Size 

In the calculations carried out thus far, we have assumed that the orifice 
was very small compared to the dimensions of the baffle and therefore the 
dimensions of the wake. For orifice sizes usually employed, this assumption 
is reasonably valid. However, there may be a variety of reasons for desiring 
a larger orifice on any particular flight. For example, if the spin rate is high 
and the chamber is large, a larger orifice area may be necessary to maintain 
steady state flow into and out of the chamber. For this reason we have calcu- 
lated the effect of the orifice size upon the : essure modulation. 

We divided the orifice into seven strips running parallel to the occulting 
baffle, as shown in figure 9. We calculated P./P Q for each strip and found the 
average P /'P by multiplying by an area factor. 


v p o = p o>j A j • 


( 22 ) 


MOTION 


♦ 


baffle 


ORIFICE 



Figure 9-The effect of the baffle occulting 
a finite size baffle. 
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For each strip, the integration limit 9 is given by 


where 


0 


sin 


RiCr 

j 


[(R ± C. r ) 2 + (i cos qp) 2 ] 12 


(23) 


R + C. r when 0 < cp < 180° 

R - C. r 180 < qp < 360° 


C. -constant determined by geometric relationship 


r -orifice radius- 


To illustrate the effect of finite orifice sizes, figures 10 and 11 show N 2 
pressure curves for R/£ = 0.05 and 0.1, respectively, and for orifice radii of 
0, 0.2 and 0.5 inches. From equation 23 it is clear that the actual sizes of R 
and £ are important for this effect, rather than just the value of R/£. We as- 
sumed a baffle distance of 10 inches and a half width of one inch and one half 
inch, respectively. 

From the figures it is apparent that the orifice radius does not significantly 
affect the depth of the pressure well until the radius approaches 1/2 inch. This 
can be understood qualitatively from figure 9. by considering the projection of 
the baffle as it occults the orifice. The initial occultation affects little of the 
total area, so the pressure reduction is not changed much. Thus the fluxes 
through the inner slices tend to dominate the shape of the pressure well. 

From this viewpoint, alternate orifice shapes might be considered. For 
example, it is apparent that only the orifice dimensions perpendicular to the 
baffle scan direction can cause important modification of the pressure charac- 
teristic. Therefore one might employ a long rectangular orifice parallel to the 
baffle to obtain a large gas conductance without reducing the depth of the 
pressure well. 
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CHAMBER PRESSURE (ARBITRARY UNITS) 



Figure 10-The effect of a finite sizeorifice (radius 0.5 inch) upon the pressure well 
for R/'t = 0.05 and employing N 2 at 1000°K. 
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CHAMBER PRESSURE (ARBITRARY UNITS) 



Figure 11— The effect of finiteorifice sizes (radius = 0.2 and 0.5 inches) upon the 
pressure well for R/'t = 0.1 and employing N2 at 1000° K. 


APPENDIX C 5 

The Neutral Wake Experiment on a Non-spinning Spacecraft i 

In this paper we have limited our discussion to a baffle mounted rigidly in 
front of an orificed-chamber on a spinning spacecraft. We then rely on the spin 
motion to scan the wake of the baffle across the orifice. Essentially the same 
effect can be obtained on a fully stabilized satellite by passing the baffle in front 
of the orifice in some suitable manner. Depending upon how the baffle move- 
ment is achieved, the theory developed here might require little, if any, 
modification. 

One arrangement which offers promise is a propeller type arrangement in 
which one or more baffles are mounted on a spinning shaft such as shown in 
figure 12. The advantages of this system are that (1) multiple R /i ratios could 
be employed to optimize the pressure characteristics for different masses or 
temperatures and (2) optimum occultation rates can be achieved by selection, or 
variation, of the drive shaft speed. It would, of course, be necessary to arrange 
the drive shaft and baffle supports so as to minimize backscattering of particles 
into the orifice. % 
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Figure 12-Neutral wake experimental arrangement on a non-spinning spacecraft. Multiple 
baffles could be used to optimize resolution for various masses and temperatures. 
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